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Abstract 

The dynamics of charge qubit in double quantum dot coupled to phonons is investigated theo- 
retically. The static bias is considered. By means of the perturbation approach based on unitary 
transformations, the dynamical tunneling current is obtained explicitly. The biased system displays 
broken symmetry and a significantly larger coherence-incoherence transition critical point a c . We 
also analyzed the decoherence induced by piezoelectric coupling phonons in detail. The results 
show that reducing the coupling between system and bath make coherence frequency increase and 
coherence time prolong. To maintain quantum coherence, applying static bias also is a good means. 
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1. Introduction 

Since the discovery that quantum algorithms can solve certain computational problems 
much more efficiently than classical onesjlj], attention has been devoted to the physical im- 
plementation of quantum computation. Nanofabrication technology now allow us to design 
artificial atoms (quantum dots) and molecules (coupled quantum dots), in which atomic 

nnrt 

(molecular) -like electronic states can be controlled with external gate voltages[2|, y, |4J. As 
the system always interacts with its environment, quantum decoherence in the system usu- 
ally is the most serious obstacle produce efficient qU an tUm c ir c„ itS fl,flQ . For the reason, 
a detailed understanding of quantum decoherence in open system and implementing suffi- 
ciently high number of coherent manipulation within the characteristic coherence time of 
qubits are crucial for future actual implementation of quantum nanostructures to quantum 
information technology. 

To build a quantum computer, the first step is the realization of the basic device units 
for quantum information processing called quantum bit (qubit). Within the last decade, 
various schemes have been proposed and many of them have even been realized, such as 
superconducting flux qubit ja, a, ID, 11] and solid charge qubit [si, |j, 12], including single 



and double dot qubit. Among them, the gate voltage controlled semiconductor charge qubit 
has the potential advantages of being arbitrarily scalable to large system and compatible 



with the current microelectronics technology. Recently, Hayashi et al.[13| have successfully 
realized coherent manipulation of electronic state in double-dot system implemented in a 
GaAs/AlGaAs heterostructure containing a two dimensional electron gas. The damped 
oscillation of population inversion is observed in the time domain and the dependent of 
decoherence rate T 2 _1 on the energy offset e is presented. In another similar experiment [l^] , 
the base material used of the charge qubit is an industry-standard silicon-on-insulator wafer 
with aphosphorous- doped active region and all operations (initialization, manipulation, and 
measurement) are achieved by capacitively coupled elements. The change in gate voltage 
V 3 3 can principally control static bias e. From above, we can see it is necessary to investigate 
the effect of the static bias on the tunneling current of the charge qubit. 

In this work, we study coherence dynamics of double QDs charge qubit using spin boson 
model with static bias, which is investigated without applying the Markov approximation 
to the electron-phonon interaction. A simple explicit expression of population inversion 
or tunneling current is presented through perturbation treatment based on unitary trans- 
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formations. The coherence- incoherence transition critical point a c versus bias is gained. 
Furthermore, the piezoelectric potential phonons induced decoherence are investigated in 
detail and possible means for maintaining quantum coherence are expressed. 

The paper is organized as follows: in Sec. 2 we introduce the model Hamiltonian for 
static bias spin-boson model and solve it in terms of a perturbation treatment based on 
unitary transformations. We analyze the result and provide proposition for how to maintain 
tunneling current in Sec. 3. Finally, the conclusion is given in Sec. 4. 

2. The model and theory 

The model we study is the spin-boson Hamiltonian for static bias jig where the two state 
system is linearly coupled to a continuum of harmonic oscillators, i.e: 

H = H S + H b + H i (1) 

where H s is the Hamiltonian of the system, of the bath and Hi of their interaction that 
is responsible for decoherence. Here 

ft = -^ + ^ (2) 
H b = ^2 ^kb^h (3) 

k 

H^-J^gM + h)** (4) 

k 

with (jj being pauli spin matrices. A describes the tunneling coupling between the two 
states while e is the energy offset between the uncoupled charge states. b^ib^) and u k are 
the creation (annihilation) operator and energy of the phonons with the wave vector k. g k 
describe the electron-phonon coupling strength. In this work we consider the static bias 
case with temperature T=0. The effect of the phonon bath are fully described by a spectral 
density: 

k 

In order that the Hamiltonian is diagonalized in a z direction, we make a displacement to 
all boson modes, 

bk = a k - J^cr (6) 
Zulu 
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where o"o is a constant and will be determined later. Then we apply a canonical transfor- 
mation, H' = exp(s)H exp(— s) with the generator 16|, [17] 



S = T,ij- k ( b t-h)(cr z -a ). (7) 

k 

Here we introduce in S a A; -dependent function The aim of the transformation is to take 
into the correlation between the two states and the bath. Thus we get the Hamiltonian H' 
and decompose it into 

H' =H' + H' 1 + H' 2 (8) 



where 

H' = -\n^a x + ^ + J2^Kh - £ |U(2 - 6) + E £^o(l - (9) 



#i = \ E^ 1 - + 6 *)('« - *o - <2 £*» E - ( 10 ) 



with 



H' 2 = -^(cosh(E-6(^-M-^) 

k 

-^(sinh(E " 6*)) - V E -Ubt - 6,)) (11) 

fc k 



^ = exp(-E#7^) ( 12 ) 

and 

e'=e-ra ,r = E-(l-6) 2 (13) 

Obviously, H can be solved exactly because in which the spin and bosons are decoupled. 
Then we can diagonalized H by a unitary matrix U, 



. u v . 

U= | | , (14) 

v —u 
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with 

1 i 1 i 

u = —=(1 - sin6)^,v = —=(1 + sinfl)^ (15) 
v 2 v 2 

and sintf = e'/W with W = (e' 2 + r] 2 A 2 ) 1/2 . 

H' is transformed as follows (to the second order of (?*,): 

H" = U + H'U = H'o + H{ + Hi (16) 



H'o = -\wa z + £ ^Kh ~ E ^*(2 - 6) + E ^^(1 - 6) 2 (17) 



E ^ - ^)(^ + 6 *) + <! £*» E ~^ h i - b ^ ( 18 ) 

U'l and i?2 are treated as perturbation and they should be as small as possible. For this 
purpose <T and are determined in such a way 

ao = -±.^ k= Uk (19) 

that 

ff l' = vE ^ IKi-V* + i(7 y) + fo fcK - ia y)] ( 20 ) 
^ ^ u k 

and i?i Ig 1 ) = 0. This is the key point in our approach. We should remark that 

£ ' =£( l + £);sin S = £(i±£. < 21 > 

So 6 is in the range ofO < 9 < n/2 and the sign of sin# is identical with static bias. If 
= 0, the model correspond to zero bias case. 

We denote the ground state of Hq as \g) = |si)|{0fc}), and the lowest excited states as 
l s 2)|{0fe}), |si)|{lfc}) where |si), \s 2 ) are eigenstates of a z (a z \si) = |si), a z \s 2 ) = — |s 2 )), 
|{nfc}) means that there are n k phonons for mode k. Thus, we can diagonalize H" as: 

H" = -\w\g){g\ + ^2E\E)(E\ + terms with higher excited states (22) 

E 
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the experiments in Ref.13,14 are performed at lattice temperature below or about 20 mK. 
At such a low temperature, the multiphoton process is weak enough to be negligible. So we 
can diagonalize through the following transformation 16j 



h>|{o fc }> = ;>>(£)|£> 

E 

l*i>|{u}> = X>(tf)|tf> 



(23) 
(24) 



where 



\E)=x(E)\s 2 )\{O k }} + Y,yk(E)\ Sl )\{l k }) 



x(E) = (l + J2 



V k 2 



(E + \w-u k y- 



(25) 



(26) 



Vk{E) 



V, 



-x{E) 



(27) 



E + \W-uj k 

with V k = rjAg k ^ k /uj k . E are the diagnalized excitation energy and they are solutions of the 
equation 



V k 2 



k ~ • 2 

So Hamiltonian can approximately be describes as: 



E + \W — cu k 







(28) 



H" = ~W\g)(g\ + ^E\E)(E\ 



(29) 



The population inversion can be defined as p(t) = (ip(t) \a z \i/j(t)), where \ip(t)) is the 
total wave function in the Schrodinger picture. Since the initialization of the charge qubit 
is used to in the state \L) , it is reasonable to choose initial state \ip(0)) = e~ s \L) \0 k ) . Then 
we can obtain 



p(t) = U(o) 



Ue iH tU+e'UU+aMU+e'Ue 



iH t 



(30) 



— sin — ip' sin 6 y* k '(E')y k (E)x(E)x*(E)e< E - E '^ 

kk'EE' 

-w 2 sin#^|x(£)| 2 \x(E')\ 2 e^-z'V 

EE' 

+uvcos6Y,ME)\ 2 e l(E+ ^ t + \x(E)\ 2 e- l{E+ ^ t ) 



(31) 
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For the different terms appearing in Eq.(31), we employ the orthogonal property 

J2yk(E)y k (E') = 5(E - E') - x{E)x{E') (32) 



k 

and the residue theorem 



I <E) | 2 e iEt (33) 



2m 



e 2 7T~ / TF2 ( 34 ) 



= e- 1 ^— / _ 35 

v 2 

Denoting the real and imaginary part of UJ _ LO ^ ±i0 + as R(a?) and =py(o;), respectively, we 
can get 

R(u,) = p y £-%- (36) 

k 

oo ; 

= foA) a p / cfc/ ? ffil- ., 2 (37) 



o 



7 H = ]T Vf 5(- - **) = ^) 2 77JT^ (38) 

where p stands for Cauchy principal value, and J{oS) is the spectral density. 

The contour integral in Eq.(35) can be proceed by calculating the residue of integrand 
and substituting in Eq.(31), we obtain 

pit) = - sin 6 + sin 6(1 + sin 6)e~ 2 ^ )t + cos 2 9 co8(u t)e^ u)t (39) 

where u>o is the solution to the equation 

uj-W - R(w) = (40) 

Thus a rather simple expression for the dynamical tunneling is obtained analytically. It 
should be noted here that for t— >• oo, our result tend towards thermodynamics equilibrium 
value [l8|, which is modulated by static bias 



Comparing with Markovian approximation, our result not only give the long time limit 
behavior, but the short time damping coherence oscillation. That must holds most promi- 
nently if semiconductor quantum dots are to be used as basic building blocks for quantum 
information processing where the operation completely relies on the presence of coherence. 

Finally, the tunneling electron population in the right dot at time t can be obtain 
n(t)=l±» = i±|W. 

Until here, our presentation is not restricted by the form of the spectral density and 
can be extended to all kinds of baths. Previous work states that, at zero temperature in 
double-dot system of GaAs material, the dominant contribution of phonons to QDs come 
from the piezoelectric coupling and the deformation potential coupling is small enough to 



be ignored. So we will use the piezoelectric coupling spectral density 
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19|: 



J{uj) = au(l - — sin— )9(u c -u) (42) 

where a is the dimensionless coupling constant, u c = s/l and u>d = s/d (s is the sound 
velocity in crystal, 1 is the dot size and d is the center-to-center distance between two dots) 
and 9(x) is the usual step function. 
3. The result and discussion 

From Eq.(39) we have been able to obtain the population inversion of charge qubit p(t) 
as the analytical damped oscillation form with frequency Uq and damping rate j(u)). In the 
following calculations, uj c is taken as the energy unit. We choose the quantum dot size 1 as 
lOOnm (approximate size for the dot in Ref.13), i.e. u c = 32.5 /xeV (or 0.05 ps _1 ). Assume 
the distance between two dots is sufficiently large, d=667 nm, correspondingly = 0.15u; c . 
The typical value of tunneling barriers in experiment [l3] are A = 9 /xeV. Without special 
indication, A is choosen 9 /ieV. 

First we illustrate the population inversion as function of time (u c t) in fig. 1(a) with five 
different static bias: e = — 8 /xeV (dashed line), 0/ieV(solid line), 8/ieV (dotted line) for 
fixed a = 0.04. It is clearly shown that the population inversion exhibit damping oscillation 
and symmetry broken. The asymptotic value of t he long time limit p(oo) is determined by 
the static bias. If bias is zero, p(oo) = 0. 

Fig. 1(b) presents the long time limit of population inversion or the thermodynamics 
equilibrium value p(oo) = —e'/W as functions of bias e. Comparing with the results of 
NIBA 



181 ] . which is expressed as p(oo) = —ej \Je 2 + A^, on the scale of the figures, the two 



curves are overlapped if the identical Ohmic spectral density J{uS) = 2au)6{uj c — uj) is used. 
While the system is in the piezoelectric potential bath, p(oo) deviate a little. 

Fig. 1(c) sketches the coherence- incoherence critical point a c versus e/u c for fixed A = 9 
\ieV in Ohmic bath. We see that bias increases a c significantly and induce a transition 
from strongly damped incoherent to coherent oscillation behavior. The similar result is also 



obtained by Klaus 



15| from Monte Carlo simulations. 



The quantum coherence depends on two factors: oscillation frequency ujq and decoherence 
rate T 2 _1 . It is already noticed that long coherence time is favorable for quantum manipula- 
tion. However low oscillation frequency means that the number of quantum operation that 
can be achieved within the coherence time is very limit [3]. So we must give attention to two 
aspects coherence time and oscillation frequency. In what follows, we analyse the effect of 
the static bias e, the tunneling energy A and e-p coupling constant a on the decoherence rate 
7 and the damping oscillation frequency u and elucidate the advantage and disadvantage 
for changing characteristic energy A and e or decreasing e-p coupling to remain quantum 
coherence. 

In contrast with experiment, we present the decoherence rate 7 (T^ 1 ) as function of 
energy offset e[Fig. 2]. The electron-phonon coupling constant a = 0.02 — 0.07 is used to 
explain the inelastic current in GaAs/AlGaAs heterostructure DQD samples 19j. In Fig.2(a) 
we shows the damping rate 7 as a function of static bias e with three different e-p coupling 
a = 0.04 (solid), a = 0.08 (dashed) a = 0.12 (dotted). These curves show clearly that 
the decoherence rate raise with increasing of the absolute value of the bias e. The vary of 
7 vs e is agree with the result of experiment. If e = /ieV and a = 0.04, decoherence 
rate 7 approximately is 0.2 ns _1 . Increasing the e-p coupling to a = 0.12 the damping rate 
contributing from piezoelectric phonons grow up to 0.50 ns" 1 , occupy more than half of the 
experimental result. It proves that the coupling to the piezoelectric potential phonons is 
one of the main decoherence mechanisms in such a double-dot system. 

The damping oscillation frequency ujq vs bias e is shows in Fig.2(b), and the parameters 
are same as those in Fig. 2 (a). The qubit oscillation frequency can be changed continually 
by bias voltage and behave nonlinear response versus e, which explain the experimental 
result in Ref.14. Setting large static bias can efficiently improve oscillation frequency. When 
a = 0.04 and e = /ieV, the frequency approximately is 2.1 GHz, which is extraordinarily 
agree with the fit result of the experimental data in Ref.15. The decoherence rate and the 
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oscillation frequency reveal symmetry about bias. 

Fig.3(a) presents the decoherence rate as function of tunneling coupling A at three dif- 
ferent static bias e = /ieV (solid), 3 /ieV (dashed), 6 /ieV (dotted), fixing a as 0.04. 
Enlarging tunneling coupling A between two dots make damping rate 7 gradually increase 
as shown in Fig.3(a). So quantum coherence can be remained by reducing tunneling barrier. 
In experiment, the barrier tunneling can be principally defined by the materials and the 
geometrical constriction between the dots in the fabrication of the device. But it is possible 
to modify the tunneling barrier by gate voltages, just as shown in Ref.13 Fig.2(d). From 
Fig. 3 (a), if A in the order of magnitude of /ieV, the decoherence rate is probably in the 
range of 0.2 ~ 1 ns _1 , which is agree well with the Ref.13. A plot of oscillation frequency u 
as function A is shown in Fig.3(b). At zero bias e = 0, the curve behavior is approximately 
line. Increasing static bias, u appear nonlinear response. Comparing two similar experiment 
in Ref.15 and Ref.16, the characterized energies E* 40.5 neV corresponding to angular fre- 
quency uj = 62 MHz in the silicon two-level quantum system are 1/37 smaller than quantum 
level spacing 1.5 fieV corresponding to oscillation frequency uj = 2.3 GHz in the reports 
for semiconductor double quantum dot charge qubit. From the calculation result, we can 
deviate that the reducing of energy split make damping rate and oscillation frequency fall 
same order of magnitude, correspondingly, in the piezoelectric coupling spectral density, so 
the lack of piezoelectric coupling in the silicon QD might result in the long coherence time 
two orders of magnitude longer than reports for semiconductor QD. We find that the bias 
can effectively enhance oscillation frequency but have relative little lifting for damping rate 
in the case of weak tunneling couple. Therefore adjust static bias larger is another choice 
to preserve coherence and that is a suggestion to overcome the obstacle in Ref.16. 

Fig.4 presents the decoherence rates 7 and oscillation frequency uj as functions of dimen- 
sionless coupling constant a at three different bias e = /ieV (solid), 3 /ieV (dashed), 6 fieV 
(dotted). As shown in Fig. 4 (a), for determined bias, it is an almost linear relation between 
damping rate 7 and dimensionless coupling constant a. Fig.4 (b) displays the oscillation 
frequency uj versus coupling constat a. These curves show clearly that when the coupling 
constant a become weaker three or two orders of magnitude from 0.04, the oscillation fre- 
quency uj have linearly increased, albeit very relaxedly. Applied static bias, the decoherence 
rate and oscillation frequency always are enhanced. As a result, to keep the coherence of 
charge qubit in long time, one good way is to minish the coupling with environment, either 
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by find good material with small e-p coupling, or by modifying the design of the DQD 
structure, or by better gate tuning. 

Anyway, to realize quantum computation, the important direction are to reduce the 
coupling to environment by all means, but to ensure effective interdot coupling and static 
bias. 

4. Conclusion 

We studied the charge qubit dynamics with static bias by a perturbation treatment 
based on two times unitary transformations. Our approach applies to all forms of spectral 
density. In Ohmic bath, the result shows that a c clearly increases and population inversion 
breaks symmetry in biased system. Analyzing the piezoelectric coupling phonon induced 
decoherence, we find that, weak coupling of the charge qubit to the environment carries 
out large coherence oscillation frequency and long coherence time. When tunneling coupling 
between two QDs is fixed, adjust gate voltage to enlarge static bias make oscillation frequency 
observably increase while damping rate unnoticeably increase. That is another better choice 
to maintain quantum coherence. Finally we hope that our predictions can offer advice for 
experimenter and be testified by experiment in the near future. 
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FIGURES 



Fig.l.(a) The population difference as a function of u c t for different static bias —8 fieV 
(dashed line),0 p,eV (solid line), 8 /ieV (dotted line), (b) The long time limit p(oo) 
(solid) as a function of static bias e in piezoelectric potential bath. p{po) (dotted) in 
Ohmic bath and p(oo) (dashed) in Ohmic bath for NIBA. The coupling constant is 
fixed as a = 0.04. (c) The coherence-incoherence critical point a c versus static bias 
e/oo c with A = 9 /ieV. 

Fig. 2. The damping rate 7 and oscillation frequency u versus e in (a) and (b), respectively 
with A = 9 /ieV, ujd = 0.15cj c for different a = 0.04 (solid line), 0.08 (dashed line), 
0.12 (dotted line). 

Fig. 3 The damping rate 7 and oscillation frequency u versus A in (a) and (b), respectively 
with a = 0.04, oo d = 0.15cj c for different static bias e = /ieV (solid line), 3 /xeV 
(dashed line), 6 /ieV (dotted line). 

Fig. 4. The damping rate 7 and oscillation frequency uj versus a in (a) and (b), respectively 
with A = 9 /xeV, LOd = 0.15c<j c for different static bias e = /ieV (solid line), 3 /xeV 
(dashed line), 6 /xeV (dotted line). 
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